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We present a study of entropy transport in Bi2Se3 at low temperatures and high magnetic fields. 
In the zero-temperature limit, the magnitude of the Seebeck coefficient quantitatively tracks the 
Fermi temperature of the 3D Fermi surface at T-point as the carrier concentration changes by two 
orders of mag nitude (10 17 to 10 19 cm -3 ). In high magnetic fields, the Nernst response displays giant 
quantum oscillations indicating that this feature is not exclusive to compensated semi-metals. The 
analysis of quantum oscillations quantifies the magnitude of the Zeeman energy and points to a 
significant field-dependence of the Fermi energy across the quantum limit. 



The B12X3 family (X= Se, Te) is attracting tremen- 
dous attention as a Topological Insulator (TI). Recently, 
the existence of this new class of bulk insulators was pre- 
dicted and confirmed [T]. In a TI, the bulk energy gap is 
traversed by spin polarized surface states. Therefore, the 
electrical conduction is expected to occur only at the sur- 
face. In practice, however, these materials are often low- 
density bulk metals. Interestingly, many of the TI of the 
first and the second generation are well-known thermo- 
electric materials [2] and present a sizeable thermoelectric 
figure of merit. This quantity, ZT = ~ (here S is the 
Seebeck coefficient, n is thermal conductivity and p is re- 
sistivity) characterizes the thermoelectric efficiency of a 
material. To this date, the largest thermoelectric figure 
of merit in a bulk material at room temperature has been 
reported in Bi 2 Te 3 (ZT=0.8 at T=300K [IS]). 

In spite of the fundamental and technological inter- 
est in Bi2Se3 and Bi2Te3, their thermoelectric response 
has never been investigated at temperatures low enough 
to distinguish between the two competing (semiconduct- 
ing vs. metallic) ground states. In this Letter, we re- 
port on measurements of Seebeck and Nernst effect in 
n-type Bi 2 Se 3 with a bulk carrier concentration vary- 
ing from 10 19 cnr~ 3 to 10 17 cm -3 at low temperature and 
in a magnetic field as strong as 32 T. We find that the 
low-temperature magnitude of the Seebeck coefficient is 
set by the Fermi temperature, which can be directly 
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extracted from the measured properties of the metallic 
Fermi surface. The transverse thermoelectric (Nernst) 
response displays large quantum oscillations. Such oscil- 
lations were previously reported in low-density compen- 
sated semi-metals, bismuth and graphite. Their obser- 
vation in a single-band uncompensated system indicates 
that they are generic to any low-carrier system pushed to 
the quantum limit by a sufficiently strong magnetic field. 
Analysis of quantum-oscillations allows us to: i) docu- 
ment a continuous field-induced shift in chemical poten- 
tial affecting the periodicity of the oscillations; ii) quan- 
tify the magnitude of Zeeman splitting in this system, 
which is expected to host a large spin-orbit coupling. 
These two points are crucial in discussing the Landau 
spectrum of the system at high magnetic field. 

Measurements of longitudinal (S=^7p) and the trans- 
verse (N= \j\ v t\ ) thermoelectric response were performed 
on a standard one-heater-two-thermometers setup. For 
all samples, the electric current and the thermal gradient 
were applied in the plane of the quintuple layers and the 
magnetic field was oriented along the trigonal direction 
(perpendicular to the layers). The samples used in this 
study are similar to those previously studied by longitu- 
dinal and transverse magnetoresistance experiments [5] 
and described there. 

In Fig.l, we compare the electrical transport and the 
entropy transport in two samples Ai and Bi with typical 
bulk concentrations of 10 19 cm~ 3 to 10 17 cm~ 3 . For both 
samples, p xx and p yx display quantum oscillations on top 
of an almost linear monotonic base as previously reported 
and discussed |5 8 . In presence of a field of about 10 T, 
we find for the two concentrations that p X x ~ Pyx 

whereas 

S >> N (the typical ratios are 100 for Ai and 10 for Bi). 
As in the case of the electrical transport, S and N display 
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FIG. 1: Electrical and entropy transport measurements of 
Bi2Se3 for the samples Ai (n(Ai) «lel9 cm -3 ) and Bi (n(Bi) 
~ lel7 cm -3 ) a) p xx (in blue) and p yx (in red) of Ai as a 
function of the magnetic field for T=1.9 Kb) S (in blue) and 
N (in red) of Ai as a function of the magnetic field for T=2.6 
K (the insert show a zoom of S and N between 8 and 12T) 
c) pxx (in blue) and p yx (in red) of Bi as a function of the 
magnetic field for T=350mK d) S (in blue) and N (in red) of 
Bi as a function of the magnetic field for T=3.1K 



quantum oscillations with a phase shift of \ . The oscil- 
lations are particularly pronounced in N (in the case of 
Ai they dominate the signal). The surprising sensitivity 
of the Nernst effect to quantum oscillations in Bi2Se3 is 
indeed reminiscent of those reported in bismuth [pj and 
graphite [ID] . The origin of these giant quantum oscil- 
lations is the subject of on-going theoretical researches. 
[PTl - fT3] . These two elemental semi-metals are, however, 
compensated systems, whereas the (bulk) Fermi surface 
of Bi2Se3, as we will see below is a single band at the 
T-point. Therefore, the observation reported here estab- 
lishes that this effect is not exclusive to compensated sys- 
tems. We note that there are several reports on Nernst 
quantum oscillations with a large amplitude and a small 
frequency in metals like zinc |14] and aluminum [15] or in 
doped semiconductors such as n-type InAs[16] and iron- 
doped HgSe [T7]. 

Next, we focus our attention on the temperature de- 
pendence of the quantum oscillations in the thermoelec- 
tric properties, which have been so far poorly studied. In 
Figl2l ^ = is plotted as a function of B _1 for various 
temperature. In both cases, the signal is periodic in B _1 . 
The difference of the period between the two samples sim- 
ply reflects the difference in the carrier concentrations. In 
both cases, the oscillations disappear at a typical temper- 
ature of 20 K suggesting rather similar low effective mass 
for both concentrations. Following the standard Lifshitz- 
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FIG. 2: a) and b) £ = ^ as a function of B -1 for the 
sample Ai and B2 (sample different of Bi with a similar bulk 
concentration) for the various temperatures explored c) Tem- 
perature dependence of the amplitude of the oscillations of 
the sample Ai (in red) and B2 (in blue). The line correspond 
to a fit using the Lifshitz-Kosevich formula (see the text ) 



Kosevich theory [18] . we fit the temperature dependence 
of the oscillating part of ^ to : RT = sin ^ x ^ where X 

= 14.69m- T ^ The data fit shown j n Fig|2]c). 

We found respectively m*(Ai)/m = 0.20 ± 0.03 and 
m* (B2) I 'mo = 0.18 ± 0.03 in good agreement with the 
cyclotron mass deduced by Shubnikov-de Haas oscilla- 
tions 0EH5J. 

The magnitude of the thermoelectric response is dra- 
matically affected by the change in the carrier concentra- 
tion. Fig|3]a) and b), present |! and ^ as a function of 
temperature for the four samples studied. As the dop- 
ing passes from 10 19 cm -3 to 10 17 cm -3 , ^ increases by 
one order of magnitude and ^ increases by two orders of 
magnitude. 

In the Fermi liquid picture for a one-band system, 
the diffusive Seebeck is expected to be T-linear in the 
zero-temperature limit with a magnitude proportional to 
1/T F : 
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Similarly the solution of the Boltzmann equation for 

the Nernst response leads to : y=-x3n^"fe 
can be simplified and rewritten as [19] : 
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Empirically, these simple equations give a rough ac- 
count of the magnitude of transport coefficients across 
several orders of magnitude [TH1 [50] . 
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TABLE I: Properties of Bi 2 Se 3 
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Bi2Se3 provides a particularly compelling opportunity 
to check the robustness of Eq. [I] and Eq. [2] in a wide 
window of carrier concentration with a barely changing 
Fermi surface topology. The mobility and the Fermi en- 
ergy can both be extracted from the data on quantum 
oscillations. Tab|l] lists the values of ^ and ^ obtained 
at the lowest temperature measured ( 2K for Ai and B 2 
and 150 mK in the case of A 2 and Bi) respectively. Val- 
ues for the Fermi temperature and the mobility extracted 
from experiment were plugged in Eq{l] and Eq[2] The 
expected values for the four samples are reported in the 
two last lines of Tabjlj The overall agreement between 
the measurement and the expected values is remarkable. 
This is particularly the case for the two samples studied 
at the lowest temperature (A 2 and Bi). 

When the carrier concentration becomes as low as 
10 17 cm -3 , the so called quantum limit can be attained 
by a magnetic field of 15 T. In this limit all carriers are 
confined to their lowest Landau level. One essential pa- 
rameter in this regime is the magnitude of the g-factor. 
In the case of bismuth, for example, when the field is 
oriented along the trigonal direction, the g-factor of the 
hole pocket is as large as 62. This corresponds to a ra- 
tio of the Zeeman energy (Ez) to the cyclotron energy 
(Hoj c ) labeled M = |f slightly larger than 2 [2"Tl 
As a consequence, the chemical potential starts to move 
significantly with increasing magnetic field, dramatically 
affecting the Landau level spectrum [22] . 

In the case of Bi2Se3, the data of sample B 2 down to 
T=0.38mK and up to B=17T does not reveal any Zee- 
man splitting of the peaks like the one found in bismuth 
[51 [21] • This suggests that the Zeeman energy and the 
cyclotron energy are commensurate (i.e that M is closed 
to an integer in the limit of the peak width) as origi- 
nally suggested by Kolher et al. [6]. In order to differ- 
entiate between M=0,l,2,3,..., we used a model similar 
to the one used by Kolher et a/.[B] that determines the 
Landau level spectrum and the field dependence of the 
Fermi energy. The calculation of E^ is performed with 
the assumption that the carrier concentration (noted n) 
is independent of the magnetic field (see part B of the 

S. I. for more details): n = J_ F D(e)de where D(e) is 
the density of state, equal to the sum of the density of 
state of each Landau level (noted LL) depending on their 
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FIG. 3: a) | as a function of the temperature for A 2 (same 
batch as Ai with a slightly lower concentration) and Bi b) & 
at B=1T as a function of the temperature for Ai, A 2 ,Bi and 
B 2 



position relative to the Fermi energy (see part F of the 
S. I. for more details). One complication arises in com- 
paring the experimental Nernst peaks and the theoreti- 
cal Landau spectrum. According to a recent theory by 
Bergman and Oganesyan[TT], it is the off-diagonal ther- 
moelectric conductivity a xy = Pm p 2 N +p2 V times magnetic 
field peaks when a Landau level coincides with the chem- 
ical potential. In the case of bismuth and graphite, since 
Pxx >> Pxy and S <N, this quantity peaks when N 
peaks. However, in the case of Bi 2 Se3 p xx N and p xy S 
become comparable in amplitude and therefore a xy x B 
and N do not peak exactly at the same magnetic field. 
However, the difference between the peak positions in 
a xy x B and in N is small compare to the distance be- 
tween two peaks in a xy or in N. As discussed in part E 
of the S.I. section, our conclusion is not affected by this 
choice. 

The last peak resolved in the Nernst effect in sample 
Bi occurs at 10. 4T. We have adjusted the carrier concen- 
tration in order to find the best agreement between the 
position of this peak and the crossing points between a 
Landau level and the Fermi energy for different possible 
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FIG. 4: In red line is the Nernst effect, N, as a function of 
B _1 at T=3K for the sample B2 on a) and b) and as a function 
of B at T=5K for the sample B3 on c) and d). The evolution 
of the Fermi energy (Ef) is plotted in black points. Each 
time that the chemical potential is crossing a Landau level 
(colors lines), there is a kink due to the carrier variation. The 
best parameter to describe the peaks position as measured in 
N is found to be for M=l (with a bulk carrier concentration 
n=1.28 10 17 cm -3 ) on a) and M=2 ( with a bulk carrier con- 
centration n=1.08 10 17 cm -3 ) on b) . c) and d) same as a) 
and b) with the sample Bs( same as B2 few months after the 
a) experiment). 



values of M. For M=0 and M=3 (and above), we did not 
find a set of parameters matching the position of all the 
peaks. However, as shown in Figj4] a) and b), a good 
agreement is found for M=l and M=2. Note that the 
agreement is slightly better for M=2 than for M=l. In 
the case of M=l, the LL n+l~ and n + are degenerate 
except for the last ones: CP (in orange on Figfila)). In 
the case of M=2, the LL n+2~ and n + are degenerate 
except for the two last ones 1~ (in green on FigHb)) and 
CP (in orange on Figjijb)). Remarkably, in both cases 
above 10. 4T, all bulk carriers are spin polarized. 

Sample B3 was studied up to 32T. In this case, the 
last resolved peak in the Nernst effect occurs at 11. 4T 
because of a slightly higher carrier concentration. In- 
terestingly, between 17T to 32T, an unexpected increase 
of the Nernst effect was observed. At first glance, one 
could explain this increase by the presence of additional 



Landau levels which are expected in the case of M>2. 
For instance for M=3, the last observed peak could be 
attributed to the LL 2~. However, this scenario does 
not reproduce the low field peak positions. Moreover, it 
would yield a chemical potential crossing for the l - LL 
at 19T, in contrast with our experimental data. Thus 
we conclude that the two only possible values for M are 
M=l and M=2, implying respectively a g-factor of 14.5 
and 29 (for an in-plane cyclotronic mass of 0.14mo). 

The increase in the Nernst effect observed at high 
field calls for an alternative scenario. As reported on 
FigjZ] b), once the l - LL crosses the chemical poten- 
tial, the enhanced shift in chemical potential pushes it 
closer towards the CP LL. In the the low field regime 
the Nernst response increases as the distance between 
the highest filled Landau level and the chemical poten- 
tial becomes shorter and attains a maximum when the 
two levels coincide. Here if the chemical potential ap- 
proaches the CP LL without crossing it in order to keep 
charge neutrality, one expects it to peak without attain- 
ing a peak. This scenario can qualitatively explain the 
experimentally-observed increase in the Nernst effect. 

In conclusion, the overall thermoelectric and thermo- 
magnetic response in Bi2Se3 can be quantitatively under- 
stood in the context of a single 3D band in presence of 
large Zeeman splitting. Our investigation reveals a large 
variation of thermoelectric and thermomagnetic coeffi- 
cients with the bulk carrier concentration which can be 
quantitatively understood as a result of a change in the 
mobility and the bulk Fermi energy. We resolved large 
quantum oscillations in the Nernst coefficient suggesting 
that such an effect is a common feature to both compen- 
sated and non compensated low carrier systems. For the 
lowest carrier concentrations investigated, we found that 
the Zeeman energy is either equal to one or two times the 
cyclotron energy. In both strong variation of the 

chemical potential above the quantum limit is expected, 
which can naturally explain the observed increase in the 
Nernst response. This field-induced shift in chemical po- 
tential should be taken into account in discussing the 
contribution of surface states to high-field transport in 
this system. 
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Supplementary Information 



A. Sample growth 



Single crystals of Bi2Se3 were prepared by melting high purity bismuth (6N) and selenium (5N) in sealed quartz 
ampoules as previously reported [SI]. The carrier density was controlled by varying the initial Bi:Se ratio. 



B. Hall effect and carrier concentration 
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FIG. 5: a) Field dependance p yx for Bi (T=1.5K up to 12T), B 2 (T=0.35K up to 17T) and B 3 (T=1.6K up to 30T) b) Field 
dependance p yx for Ai and A 2 at T=1.6K up to 12T 



Figj5] shows the field dependence of the Hall resistivity for five samples. Up to 32T and down to the lowest 
carrier concentration investigated, the Hall effect was found to be linear in magnetic field suggesting that the carrier 
concentration remains constant and does not vary with the field. From the slope, we can determine the carrier 
concentration, labeled n H and reported in Tabjll| These values can be directly compared with the carrier concentration 
deduced from the period of the quantum oscillations labeled n F in Tabllll The conduction-band structure in this 
range of carrier densities is approximately parabolic. The Fermi energy is then determined by the period of the 
oscillation through the Onsager relation. Assuming a mass anisotropy independent of the doping ( mi=m2=0.18mo 
and m.3=1.7mi) one can determine the carrier concentration of the samples. As seen in the table, for all samples, the 
carrier concentrations determined by the two methods are in good agreement. 
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TABLE II: Carrier concentration and electronic mobility deduced from various transport measurement. n H and n F correspond 
respectively to the carrier concentration deduced from the Hall effect and from quantum oscillations of the Nernst effect. 


Sample Ai° 


A 2 " 


Bi 6 


B 2 a 


B 3 a 


n"{cm- J ) 3xl0 19 
n F (cm- 3 ) 2.4 xlO 19 


1 x 10 iy 
1.5 x 10 19 


2.5 x 10 1Y 
3.8 x 10 17 


4 x 10 1Y 
4.8 x 10 17 


3 x 10 iY 
4.4 x 10 17 


p(Magnetoresistance) (cm 2 .V -1 .s) 500 
hIhqII) (cm 2 .V _1 .s) 670 


800 
1200 


4500 
4740 


5000 
5330 


^4500 
r;4000 



"measured down to 2K 
'measured down to 0.15K 

C. Electronic mean free path and impurity scattering 

In the case of Bi2Se3 it is believed that Se vacancies are the charge dopants responsible for the n-type nature of 
the system. Here we propose to make a simple discussion to know if the vacancies are the main source of scattering 
of the electrons at low temperature. In other words we propose to compare the electronic mean free path and the 
distance between vacancies (d v ). For the two extreme concentration that we studied, we can determine the electronic 
mean free path l e using the knowledge of the mobility and the Fermi velocity : l e =vpT= hkp ^ . 

Each Se atoms can contribute for two electrons. For a given bulk carrier concentration n, the number of lacuna per 
unit cell, Ni is given by : N+= where V the unit volume. If we assume a homogeneous distribution of the lacuna, 



the typical distance between vacancies sites can be estimated : |7rd^= = ^. Tab III yields values of d v for the two 



carrier concentrations. For both samples, the mean-free-path is longer than the typical distance between impurity 
sites, suggesting that the scattering probability between a traveling electron and a point defect is much smaller than 
unity. 

TABLE III: Comparison of the electronic mean free path and the distance between vacancies sites 



Sample 


n (cm 3 ) 


le(nm) 


d v (nm) 


Ai 


3 x 10 19 


23 


2 


Bi 


2.5 x 10 17 


550 


10 



D. Dingle Temperature 

From the temperature dependence of the amplitude of the Nernst coefficient we extracted the in plane effective 
masse. By studying the held dependence, we can also extract additional information such as a Dingle temperature. In 
Fig|6j we report the so-called "Dingle plot" for four temperatures T=2.8 K, 8.3 K, 12.5 K and 15.8 K for the sample 
Ai . In the context of the Lifshitz-Kosevich formalism [S2] , we expect that : 

log(^) = -aT D m* x I (3) 

For the four temperatures, the field dependence of log(^) is linear. From Eq 3 we found a Dingle temperature 
Td=10 ± IK in good agreement with Shubnikov-de Haas measurements [S1,S3,S4^S5]. 

E. 

Several theories were proposed to explain the observation of giant quantum oscillations in bismuth and graphite 
[S6,S7,S8]. Two of these theories focus on the off-diagonal thermoelectric conductivity a xy [S6,S7], another one is 
based on a semi classical description of the Nernst effect (N) [S8]. In the case of bismuth and graphite, these two 
descriptions are equivalent because p xx » p xy and S<<N. However, the equivalence between a xy * B and N does 
not hold in the case of Bi2Se3 since p xx rj p xy and S> N. We report on FigJT] a) the field dependence of the four 
quantities: a xy * B, N, Ap xx and Aa xx —cr xx (0.25K) — a xx (A.2K) (where a xx = 2 P T 2 ) as a function of the magnetic 

field. Several conclusions can be drawn from this plot: 
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First, a xy * B and the Nernst effect do not peak exactly at the same magnetic field. Fig[7]b) shows the B _1 position 
of the Landau levels vs. the index number. As seen in Figj7]a), the difference between the peak positions in a xy * B 
and in N is very small compare to the periodicity of the signals. Thus, our conclusions is not affected by the choice 
of a xy * B or the Nernst effect. Further theoretical investigations are needed to clarify this issue. 
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FIG. 7: a) Field dependence of the Nernst effect at T=3K, a xy * B(T=3K), Ap xx and Aa xx =a xx (0.25K) - a xx (A.2K) 
respectively in red, green, blue and yellow, b) Landau-level fan diagram for oscillations in N and a xy * B respectively red and 
green. The red line corresponds to a linear fit of the red points. As seen in a), the peak positions in N and a xy * B are almost 
merging 

Second, the maxima of Aa xx are concomitant with the peak positions in a xy and the Nernst effect, whereas the 
maxima in Ap xx are not. The origin of this difference is again related to the similar amplitudes of p xx and p xy above 
a few Tesla. Thus it appears that the enhancement of the conductivity generated by the crossing of a LL (Landau 
level) and the Fermi energy leads to a maximum in the Nernst response in the case of Bi2Se3, as it is the case for 
graphite and bismuth. 

In Figj7]b) (Landau level fan), the intercepts of a xy *B and the Nernst effect are very close to 0, which is reminiscent 
of the non trivial Berry phase as observed in LaRhlns [S9]. In the case of Bi2Se3, which is characterized by a parabolic 



dispersion, the large spin orbit interaction generates this peculiar behavior. The occurrence of a vanishing intercept 
points out to an odd value for the ratio of the Zeeman energy and the cyclotronic energy (labelled M). This value is 
different from that deduced from the analysis of the peak position of Ap xx by Kohler et al. [S3] (M=2). However, 
as discussed previously, it has been shown that the peak position in Ap xx differs from the Nernst response (and 
also differs from a xx ). In addition, as we approach the quantum limit, the spectrum analysis becomes more difficult 
because of the field dependence of the Fermi energy. We propose now to describe the theoretical background used to 
pin down the possible values of M. 



F. Zeeman energy in Bi2Se3 
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FIG. 8: Nernst effect (red line) of the sample B2 at T=3K on top of a calculated Landau level energies and Fermi energy 
Ef (black open circle) with different parameters: M=0,l,2 and 3 which respectively correspond to a), b), c) and d). The 
calculation was performed with a carrier concentration of n=1.08 10 17 cm -3 for a), b) and d) and n=1.2510 17 cm~ 3 for b) , a 
Dingle temperature of T_d=8K and mi=m2=0.14mo and m3=0.24mo. 

In presence of a large <?-factor (M^l) for one single band, the Landau spectrum is modified in two ways: 

i) The quantum limit is reached at a higher magnetic field. For example when M=2, the quantum limit is approxi- 
matively one period above the quantum limit attained than in the case of M=0. 

ii) The lowest LL 0~ is going down with increasing magnetic field. This pulls down the chemical potential near the 
quantum limit regime. In order to quantify these effects in the case of Bi 2 Se3, we compute the field dependence of 
the Fermi energy for various value of M. As discussed in the section [B| we assumed that the carrier concentration is 
independent of the magnetic field and is given by n = f_£ Q D(e)de where D(e) is the density of state. In presence of 
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a magnetic field and in the absence of a spin splitting, the density of state Do(e) can be written as : 

1/2 



A)(e) 



/mo 



n=0 



(e - e r 



r 2 



+ r 2 



1 if e > /iw c /2 
if e < /iu; c /2 



(4) 



=(n+5)So; c , n mQX is the highest value of n yielding a positive value value for £ — £„ and T is the broadening 
< j n presence of a spin splitting e„=(n + | ± ^-)huj c where M — 
the density of state is given by : D(e) = ^(D a (e — ^-huj c ) + D (e + 4^fiw c )). 



Where e„ 
of the Landau levels (Y=ir- 



2mo 



In this case, 
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FIG. 9: N (in red lines) of the sample B2 at T=3K on top of a calculated Landau level energies and Fermi energy Ef (black 
open circle) with M=0 for a) and M=3 for b). In a case of a) the carrier concentration has been chosen in order that the last 
resolved peak corresponds to the LL 1* (n=1.810 17 cm 3 ). For b), the carrier concentration has been chosen in order that the 
last resolved peak matches the LL 2~ (n=2.910 17 cm 3 ). The calculation was performed with a Dingle temperature T_d=8K and 
mi=ni2=0.14mo and m3=0.24mo 



Using this simple model, we can adjust the carrier concentration n and the value of M to match the peak positions 
as observed in the thermomagnetic response. Fig{8]a), b), c) and d) shows the Nernst effect (red lines) as function 
of B _1 , superimposed with calculated Landau levels and Fermi energy (black points) with various parameters: 
M=0,l,2 and 3. 

a. M—0 : On Fig{8]a), we adjust the carrier concentration to match the periodicity of the Nernst response 
(n=1.08 10 17 cm~ 3 ). In this case, the crossings of the LL and the Fermi energy correspond to minima and not 
maxima. We cannot reproduce the lowest minimum of the Nernst effect. In Figj9]a), we try to adjust the carrier 
concentration to match the last resolved peaks at 10. 4T and the LL 1 . However, with this carrier concentration 
(n=1.8 10 17 cm~ 3 ), the low-field peaks cannot be reproduced at the right positions. 

b. M=l : Figj8]b) shows the Landau level spectrum for a bulk carrier concentration of n=1.2810 17 cm~ 3 . In 
this case, there is not a complete agreement for the two last peaks, but it is rather good for the other peaks. The last 
observed peak at 10. 4T is attributed to the Landau levels + /l~ and the Landau levels n-l + are degenerate with n~. 

c. M—2 : Figj8]c) shows the Landau level spectrum for a bulk carrier concentration of n=1.0810 17 CTO -3 . We 
find a rather good agreement for all the peaks. In this case, the last observed peak at 10. 4T is attributed to the 
Landau level 1" and the Landau levels n _ are degenerate with n-2 - . This result is compatible with the conclusion of 
the early Shubnikov-de Haas analysis [S3]. 

d. M=3 (and above) : We also investigated higher values for M. On Figj8]d) (M=3), we report the Landau level 
spectrum for a carrier concentration of n=1.08 10 17 cm~ 3 . The lowest observed peaks match well the 1~ LL. This 
result is rather surprising as we would have naively expected that the field scale associated with the 1~ LL would 
continuously increase with the value of M. In fact, even if the two lowest LL 1~, - and the chemical potential are 
decreasing faster for M=3 than for M=2, the crossing point between the LL 1~ and the chemical potential occurs at 
the same field for M=2 and M=3. This conclusion is true for all the values M> 2. The situation is however different 
for the other LL, where the crossing points differ for M=2 and for M=3. In other words, for M=3 we can find a 
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carrier concentration where the 10. 4T peak matches the LL 1~, however the peak positions at lower field cannot be 
explained. On Fig[9]b), we investigate the possibility that the lowest observable peak is not associated with the 1~ 
LL but with the 2~ LL. We found a carrier concentration of n—3 10 17 cm 3 . However, the peak positions at lower fields 
cannot be explained. 

In conclusion, the M=l and M=2 hypothesis yields the best agreement between the peak positions in the Nernst 
effect and the crossing of Landau levels and the Fermi energy. For M=l and M=2 and with a cyclotronic mass of 
0.14m , we find respectively that g=14.3 and g=28.6. 
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